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■ Cosmological inflation generates a spectrum of density perturbations that can seed 

CN ■ the cosmic structures we observe today. These perturbations are usually computed as 

the result of the gravitationally-induced spontaneous creation of perturbations from 
\ an initial vacuum state. In this paper, we compute the perturbations arising from 

\ gravitationally-induced stimulated creation when perturbations are already present 

PsJ ' in the initial state. The effect of these initial perturbations is not diluted by infla- 

tion and survives to its end, and beyond. We consider a generic statistical density 

operator p describing an initial mixed state that includes probabilities for nonzero 
\ numbers of scalar perturbations to be present at early times during inflation. We 

c~| I analyze the primordial bispectrum for general configurations of the three different 

Q-r momentum vectors in its arguments. We find that the initial presence of quanta can 

significantly enhance non-gaussianities in the so-called squeezed limit. Our results 
. show that an observation of non-gaussianities in the squeezed limit can occur for 

single-field inflation when the state in the very early inflationary universe is not the 
vacuum, but instead contains early-time perturbations. Valuable information about 
the initial state can then be obtained from observations of those non-gaussianities. 

; PACS numbers: 98.80.Cq, 04.62+v, 98.70.Vc 
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O ' I- INTRODUCTION 

2 ! Quantum field theory in curved spacetimes can be considered at the present time as a 

well established framework (a recent introduction can be found in [l| and references cited 
there). In this framework, quantum theory and the effects of general relativity are brought 
together in the regime where we are confident about the validity of both theories. One 
studies, then, the behavior of quantum fields propagating in a spacetime described by a 
classical metric as in general relativity. In that way, the fundamental difficulty of formulating 
quantum field theory in the absence of a spacetime background and treating the quantum 
aspects of nonlinear gravity is avoided. However, in spite of its semiclassical character, 
very profound physical results have been obtained in quantum field theory (QFT) in curved 
spacetimes. One of the most remarkable outcomes of this framework is the phenomenon 
of spontaneous creation of quanta by curved spacetimes, as first pointed out and analyzed 
in @, |] in the cosmological context of an expanding universe. Fundamental implications 
of this phenomenon have been obtained in physical situations where spacetime curvature 
plays an important role, namely in the vicinity of black holes and in the early universe. 
In the former case, Hawking predicted the spontaneous emission of thermal radiation by 
black holes that are produced by gravitational collapse Q. In the cosmological scenario, 
when the idea of inflation was investigated [5|, it was predicted (fl that the inflationary 
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expansion of the universe is able to spontaneously create a density of perturbations and 
gravitational waves with a nearly scale-invariant spectrum. Black hole thermal radiance 
has not been detected yet, implying a low density of radiating primordial black holes in 
the observable universe. On the other hand, the high degree of precision in measurements 
of the temperature fluctuations of the CMB and the large scale structure (LSS) achieved 
during the last decade, together with the new generation of observational missions, such 
as the PLANCK satellite, open the possibility of scrutinizing in great detail the possible 
effects of the physics of the early universe. In particular, in the inflationary scenario the 
temperature fluctuations of the CMB and the LSS are directly related to the spectrum of 
density perturbations created during inflation. Therefore, the observation of the CMB and 
LSS open a fascinating window to measure direct consequences of the predictions of QFT 
in curved spacetimes, converting cosmology into a suitable scenario where this fundamental 
physics can be tested. 

The present observation of the CMB temperature inhomogeneities [3] indicates that they 
are distributed following a nearly scale-invariant spectrum. They thus are compatible with 
the predictions of inflation. The new generation of observational missions will significantly 
improve the precision of the measurements of CMB temperatures and polarizations, and of 
the LSS, in such a way that they may constitute a test of the inflationary paradigm and 
could provide enough information to favor some specific models for inflation and rule out 
others. In particular, the study of the non-gaussianities of the temperature fluctuations of 
the CMB and LSS will offer a sharp tool to obtain specific information about the dynamics 
of the field or fields driving inflation (see [§[ and references therein). 

The predictions of inflation for the spectrum of density perturbations produced at the end 
of inflation are generally computed by studying the spontaneous generation of quanta from 
a suitable vacuum state (see for example |9l-ll2|). Such a state for describing the inflaton 
quantum fluctuations, in the context of slow-roll inflation, is chosen as the natural extension 



of the Bunch-Davies vacuum |13j. This vacuum state is compatible with the symmetries of 
the spacetime in the limit where the spacetime geometry is exactly de Sitter. The election 
of the vacuum as the state describing the inflaton quantum fluctuations at the beginning of 
inflation can seem unnatural due to our ignorance of the physics in the early universe before 
inflation. However, it is consistent with the argument that the exponential expansion of the 
universe during inflation will dilute any possible quanta present in the initial state and will 
drive any arbitrary state to the Bunch-Davies vacuum. That justifies the vacuum state as 
the most natural choice. But generically, the phenomenon of spontaneous emission of quanta 
is accompanied by the corresponding stimulated or induced counterpart. As first pointed out 
in [2I, (H , the gravitationally- induced spontaneous creation of quanta in a general expanding 
universe is accompanied by the corresponding stimulated process if there are quanta already 
present in the initial state, which most generally could be a mixed state described by a 
quantum mechanical density operator, p. The objective of this paper is to study the effects 
of the stimulated creation of quanta when a general mixed state, p, describes the inflaton 
quantum fluctuations at the beginning of inflation. 

The stimulated particle creation process (as may be seen from 0, \^) is such that the 
effects of the modes that are occupied survive to the end of inflation and beyond. The 
energies and momenta of the initial occupied modes are considerably red-shifted by the 
present time. Nevertheless, the common argument that an initial density of particles present 
before inflation started would become negligible after many e-foldings of inflation is not 
correct. On the contrary, the stimulated creation process during inflation is enhanced by 
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the same quantum amplification factor as is responsible for the spontaneous particle creation 
from the vacuum that gives rise to the non-negligible amplitude of perturbations present 
after inflation. For the same reason that the effects of the quantized field fluctuations are 
observable in the spectrum of perturbations of the CMB and in the LSS, the effects of quanta 
present at the start of inflation are not dispersed away by inflation and may possibly give 
rise to observable consequences today 

In this paper, we start with a general mixed quantum state p describing a statistical 
distribution of initial inflaton perturbations parametrizing our ignorance about the physics 
before inflation. Such initial perturbations would appear to be plausible in models in which 
our universe originated through inflation from a much larger (possibly inhomogeneous and 
non isotropic) universe. Previous works have considered the possibility of having a non- 



vacuum initial state by considering a thermal state at a given temperature |14H16| or an 
exactly pure state given by a Bogoliubov transformation of the Bunch-Davies vacuum 
191 ] similar to the state first considered in 20|. Here we consider a more generic mixed state, 
described by a statistical density operator formed from individual pure states having definite 
initial numbers of perturbation quanta (i.e., particles). The initial density operator need 
not necessarily be isotropic, within the observational constraints. The condition that the 
contribution to the energy density does not significantly influence the inflating background 
geometry largely restricts the number of initial quanta. We compute the expression for 
the two-point function in momentum space (i.e., the power spectrum) for scalar curvature 
perturbations in our general mixed initial state and find constraints on the state that follow 
from the observed values of the amplitude and spectral index of the scalar power spectrum. 
We then compute the complete expression for the three-point function in momentum space 
(the bispectrum) using the canonical Lagrangian for curvature perturbations. We find that 
the three-point function is much more sensitive than the two-point function to the initial 



state, in agreement with previous analyses [15|, [18|, [19| . In particular, we obtain a significant 



enhancement of the non-gaussianties for certain configurations of the three momenta involved 
in the three point function when the number of initial quanta is equal to or greater than 1. 

However, in contrast with previous analyses, we find a significant enhancement of the 
non-gaussianities for the configuration in which two of the momenta that appear in the 
bispectrum are much bigger than the third one, for instance k\ ~ k 2 ^> k 3 (the so-called 
squeezed configuration). In this configuration, the value of the f^i parameter that char- 



acterizes the level of non-gaussianities [2l|, |22| can be enhanced by a factor k\/k% with 



respect to its value when the pure vacuum state is used. This enhancement factor can be 
as large as a few hundred for the scales accessible in observations. Our result means that 
if non-gaussianities are found to be observable in the squeezed configuration, then single- 
field inflation is not necessarily ruled out, as would follow from j23[. We show that such 
an observation could be interpreted as a consequence of an initial state for inflation that 
contains initial perturbations that have effects observable today because of the stimulated 
production of quanta during inflation. 

The rest of the paper is organized as follows. Section II is devoted to the introduction 
of the mixed initial state that we use in this paper and to the restriction that gravitational 
backreaction imposes on it. In section III, we give a pedagogical introduction to the phe- 
nomenon of stimulated creation of perturbations in an expanding universe. In section IV, 
we work out the expression for the power spectrum of scalar perturbations when the initial 
mixed state contains non- zero numbers of quanta. In section V, we present the computation 
of the three-point function and compute the level of non-gaussianities produced in different 
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configuration. We finish in section VI with some conclusions and final remarks. 

II. CHOICE OF STATE AT AN INITIAL TIME 

In the context of inflation, the density and temperature fluctuations that we observe at 
a given scale in the present universe originated at the time when the modes corresponding 
to that scale left the Hubble radius during inflation. The modes that are just coming into 
our observable universe at the present time exited the Hubble radius around 60 e-foldings 
before the end of inflation. Any shorter wavelengths observable today would have exited the 
Hubble radius at later times during inflation. Therefore, around 60 e-foldings is the earliest 
time that we can explore through observation of the non-unifomities of the present universe. 
We take that time (or a few e-foldings earlier) as the time to for setting up initial conditions 
for the quantum state of the inflaton perturbations that we want to study. At the time to, 
the energy density of the universe was some orders of magnitude below the Planck scale (see 
for example 0] for details) and the modes that we observe today in the inhomogeneities of 
the CMB and LSS were sub-planckian at that time. Therefore, we do not need to consider 
the excitation of trans-planckian modes for studying the generation of the non-uniformities 
that we observe today, and we can use QFT in curved spacetime as a plausible effective field 
theory. 

As discussed earlier, the state of the inflaton quantum fluctuations at this initial time t is 
necessary to compute the density of perturbations produced during inflation that eventually 
give rise to the temperature fluctuations of the CMB and to the LSS. In this situation, it is 
natural to proceed by writing the most general initial state for inflaton fluctuations using a 
reasonable number of variables parameterizing our ignorance about the initial state. Then, 
we can use the observational data to constrain these variables. 

We describe the state of the fluctuation field at the time to by means of a statistical 
density operator p. The possibility of states having nonzero numbers of fluctuation quanta 
in the modes at t will be included, along with the possibility of the vacuum state, in 
the operator p. Such a mixed state at to is natural within the context of inflation. In 
inflation, our observable universe is a portion of the region of space that results from the 
enormous expansion of a small patch of the initial universe. The state of such a patch would 
naturally be a mixed state because many features of the global state of the whole universe 
would not be accessible to our observable universe. As described in [24]], the state of our 
subsystem (i.e., our observable universe) arises from tracing out the full density operator of 
the large universe over "unobservables" from the point of view of the patch giving rise to 
our observable universe. That state would be a mixed state described by a density operator, 
even if the global state of the larger universe were a pure state. In such a scenario, the 
initial density operator of our universe need not be that of an equilibrium canonical or 
grand canonical ensemble. 

Motivated by such considerations, we take a general density operator p formed from 
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individual pure states having definite initial numbers 1 of perturbation quanta (i.e., particles) 

oo oo 

P=J]J]"' P(nitii), Mh), ■ ■ •) %(j2), • • •) (ni(j'i), n 2 (j 2 ), • • • I , (2.1) 

m=0 ri2=0 

where the complete orthonormal set of basis states is given by the states 

Mix), n 2 ( j2 ), ...) = (mw ■ ■ ■ r^cvrcV)" 9 ■ ■ ■ io>, (2.2) 

containing ni quanta in mode ji, n 2 quanta in mode j 2 , ■ ■ ■ , and the probabilities p sum to 
1 

00 00 

EE--PNji),n 2 (j2),..) = l. (2-3) 

rii=0 ri2=0 

In arriving at the above expression for p, the quantized scalar perturbation field, call it tp, 
has been expanded as 

p{x, t) = Y, ( A m& *) + 4/V/(*, t)) , (2.4) 

in terms of a complete orthonormal set of wave packet solutions, <Pj{x), of the perturbation 
field equation. The index j is discrete, and the annihilation and creation operators for 
quanta of the perturbation field satisfy the commutation relations 

[Ai,Al]=6 itj , [A,^-]=0, (2.5) 

where 5y is the Kronecker delta. 

The expectation value of any Hermitian observable operator B formed from p and its 
conjugate momentum is given by 

00 

(B) = Tr[pB] =J2'-- PMiO, • • •) (ni(ji), ...\B |mC/i), • • •> • (2.6) 

ni=0 

It is clear from (12. 6 \ that, unless the operator B contains terms that have an equal number 
of creation and annihilation operators, and A, , for each mode, the expectation value of 
B will vanish. 

In the Heisenberg picture, the vector states and hence p, do not evolve, but operators 
like B carry the time dependence. The reverse is true in the Schrodinger picture, so that p 
alone carries the time-dependence. In the interaction picture, the vector states are evolved 
by the interaction part of the Hamiltonian and the measured observables like B, are evolved 
by the unperturbed part of the Hamiltonian. 

The information encoded in p that will be relevant for the computation of the power 
spectrum and the bispectrum of scalar perturbations in the following sections is the average 
numbers of initial quanta present in a given mode, that is, the quantities like Tr\pNg], where 

is the operator for the number of quanta in mode with coordinate momentum k. Here, 



A mixed state formed from coherent states having indefinite numbers of initial particles would also be of 
interest, but we limit ourselves here to the former type of mixed state. 
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we have written k because the discrete set of modes we use in next sections are plane waves 
obeying periodic boundary conditions. 

However, before considering the direct observational constraints on Tr[pNr\, it is impor- 
tant to note (see [l7, 18]) that strong conditions can be imposed over these average numbers 



of particles, or quanta of the perturbation field, on theoretical grounds. On the one hand, an 
acceptable initial state must insure that physical quantities such as the energy-momentum 
tensor do not diverge when using standard methods of renormalization in curved spacetimes 
(see, for instance, [l|). That condition restricts the behavior of the function Trf/x/V^] for 
large values of k = \k\ to fall-off as ~ l/k 4+s , with 5 > 0. 

A more restrictive condition for the value of TrfpiVg] comes from the requirement that 
backreaction effects from the perturbations present at t should not modify the inflationary 
dynamics described by the non-perturbative part of the inflaton field. That implies that the 
contribution to the energy density of the inflaton perturbations at t has to be negligible 
compared to the energy density of the unperturbed part of the inflaton field, which is given 
by the inflaton potential V(4>) ~ MpH 2 . One can obtain an estimate of the restrictions 
that this condition imposes on the value of Tr [piV^] by considering the crude approximation 
TrfpA^] « N for all k/a(r ) < p, where N is a constant, and Tr[piV^] = for all k/a(to) > p, 
where p can be considered as the energy scale up to which the effective theory that we are 
considering is valid. The initial state energy density is then 

/ d 3 kkN ~ ^p 4 iV < MpH 2 . (2.7) 
Jo a W 



a\t) 



The strongest restriction on N coming from this inequality is obtained when t = £ , when the 
previous condition requires that N < MpH 2 / p 4 . For the reasonable value p ~ y/MpH ~ 
F(0) 1/4 , it follows that N can be of order one. If we take p ~ H, larger values of N are 



allowed. A more detailed analysis, including effects of renormalization, can be done [17(| (see 
also (lEI), and similar results, up to factors of order one, are obtained. Thus, renormalization 
and backreaction impose conditions on the average number of particles at to, but do not 
necessarily restrict us to small values of TrfpiVg] for all k. In sections IIVI and IVl we derive 
expressions for the power spectrum and bispectrum, respectively, of scalar perturbations 
present during inflation for an initial state described by a density operator p of the type 
described in this section. The resulting expressions for the spectrum and bispectrum can 
be directly compared with observational data to impose further restrictions on Tr[pAV|. As 
we will see, while the power spectrum imposes minor constraints, the bispectrum is a sharp 
tool to restrict the initial state for scalar perturbations. 



III. STIMULATED CREATION OF PERTURBATIONS IN AN EXPANDING 

UNIVERSE 

The effect of the stimulated creation of particles or perturbations of a scalar field was 
first derived in an expanding universe in (1, Ejj, where the effect was given for the particle 
number of a quantized scalar field, or equivalently, for the number of quanta of a scalar 
perturbation field. (As noted earlier, we can refer to quanta of the perturbation field as 
particles of that field, since they are governed by similar linear wave equations.) Here we 
focus on the dispersion of the scalar field perturbations . For readers who may not be familiar 
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with this early work, we find it useful in this section to give a pedagogical introduction to 
the phenomenon of stimulated creation of perturbations in an expanding universe. 

Consider a spatially flat Friedmann- Robertson- Walker spacetime given by the line element 

ds 2 = dt 2 - a 2 (t)(dx 2 + dy 2 + dz 2 ) , (3.1) 

where a(t) is the scale factor. In this geometry let us consider, for simplicity and for 
physical interest (next sections), a free, massless, minimally coupled scalar field <p(x,t) (the 
generalization to the massive case and the inclusion of arbitrary coupling to the Ricci scalar 
curvature is straightforward) satisfying the Klein-Gordon equation 

Op(x,t) = Q. (3.2) 

Let us proceed to quantize the field <p(x,t). In general, the statistical density operator p 
is a positive definite Hermitian operator in a Hilbert space. Then p must have a discrete set 
of eigenvalues j26[. Indeed, the set of probabilities p(ni(ji), . . .) (i.e., the eigenvalues of p) 
defined in the previous section are discrete. We find it convenient to employ spatially periodic 
boundary conditions in the quantization, such that (p(x+nL, t) = <p(x, t), where n G Z 3 . The 
use of periodic boundary conditions introduces discreteness in the values of the momenta 
k = 2nL~ 1 n. Working with discrete momenta also has advantages at the technical level, and 
in the fact that the interpretation of quantities like the number operator is straightforward 
compared to the continuum case. One can also regard the periodic boundary conditions as 
a convenient mathematical device, such that after quantities of physical interest have been 
calculated, one can take the limit L — > oo. 2 

Consider the expansion in modes of the field operator 

<p(£, t) = £ (A m (t) + A^l{t)) e rkS , (3.3) 

k 

where At and A^ are (time-independent) creation and annihilation operators associated with 
the modes (pk(t) satisfying the standard commutation relations [A%, Ag\ = [A^^Ap] = 
0. We study the spectrum of perturbations of the quantum field <p(x, t) as a function of time 
when the state is described by a density operator p. A straightforward calculation shows 
that 3 

Tr[p^ 2 (x,i)] = Tr[p^(l + tf 5 + M 5 )Mf)| a ] = E i 1 + ™\P N d) \<PkW > (3-4) 

k k 

where iV^ = AtA^ is the number operator, and we have taken into account the definition of 
p, and the fact that Tr[p] = 1. 

A general example, in which the stimulated creation of quanta of the perturbation field is 
clearly illustrated, is a spacetime for which the cosmological scale factor a(t) has an arbitrary 



2 The continuum limit can be obtained by using the asymptotic relation (27r/L) 3 ^£ — > J d 3 k, 
(L/2ir) 3 / 2 Aj: A(k) and (L/2ir) 3 / 2 ip k (t) — s- (f(k,t), where A(k) and <p(k,t) are the annihilation op- 
erator and the field mode for continuum k. 

3 The mean value of the field, Ti[ptp(x,t)], vanishes. 
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time dependence that asymptotically approaches constant values at early and late cosmic 
times t 

a{t) J aiaS ^-°°' (3.5) 
w [ a 2 as t ->■ +00 . v ' 

We assume that a(t) is sufficiently smooth and approaches constant values sufficiently fast 
that the statements we make below are well defined. Examples that allow exact solutions can 
be found in [l|, 27j and, joined smoothly to an arbitrarily long period of de Sitter inflation, in 



28J. Because at early times the geometry approaches Minkowski spacetime, it is natural to 



consider in the expansion of the field (13 .3p normalized field modes such that they approach 



plane waves at early times, ipk(t — > —00) = e~ tWlt / \J alL 3 2w\, where t and x are coordinates 
of an inertial observer located in the past, and W\ = |fc|/ai is the frequency as measured by 
this observer. The creation and annihilation operators At and then create and annihilate 
quanta representing standard particles in the early time Minkowski spacetime, and TrfpiV^] 
corresponds to the average number of these initial particles. Then, at early times, expression 
(13 .41) reduces to 

Tr[p^ 2 (f , t -+ -00)] = £^(1 + 2 Tr[pNi) , (3.6) 

k 

where we can see the contributions coming from vacuum fluctuations and from the quanta 
that are present in the initial state. 

At late times, as a consequence of the expansion of the universe, the field modes (fk{t) 
no longer will be purely positive frequency plane-waves. In general, they will be a linear 
combination of positive and negative frequency plane-waves of the form 

where W2 = \k\/ci2 is the frequency of the mode as measured by an inertial observed located 
in the late time region. The coefficients a k and fl^ are the so-called Bogoliubov coefficients 
(see [l| for additional details). They depend on the particular form of the expansion factor 
a(t) and satisfy the relation |a fc | 2 — |/3 fc | 2 = 1. 

The late-time amplitude of the modes is given by 

\Mt -> oo)| 2 = [K| 2 + Wk\ 2 + 2^(a k ^e-^)] « ^^(H' + l^ 2 ) ' ( 3 - 8 ) 

where we have neglected the oscillatory term because it averages to and is not important 
for our present discussion. Putting f !3 .81) in (13 .4ft . one obtains 4 



k 



2|/3 fc | 2 ) + 2 Tr\pN£(l + 2\l3 k \ 



(3.9) 



4 The analogous equation for the average number of particles present at late times, Tr[pNk{t — > 00)], is 
derived in Q [see second paper, Eq. (53)], and in 0] [see Ap pen dix CII ]. This stimulated particle creation 



is confirmed numerically for many e- foldings of inflation in 28 1. 
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In the previous expression, the first term in square brackets, (1 + 2|/3fc| 2 ), has two parts, 
a 1 coming from the vacuum fluctuations of the late-time Minkowski space, and a 2|/3fc| 2 
coming from the spontaneous creation of quanta that would result from the expansion of 
the universe if there were no quanta present at early-times. The 1 coming from vacuum 
fluctuations produces a divergent contribution that is removed by renormalization in the 
late-time Minkowski space. 

The second term in square brackets, 2 Tr[p7VV](l + 2\f3 k \ 2 ), is proportional to the average 

number of particles in mode k that were present in the early-time Minkowski space. The 
first part of this term is 2 Tr[pNr], which comes from the original particles that were present 
in the early-time Minkowski space. If the universe expands by a large factor, aija\ 3> 1, 
then the density of these original particles is greatly diluted by the expansion, and this term 
has a negligible effect. The second part of this term is 4 1/3^1 2 Tr[pNr], which comes from 
stimulated particle creation. The presence of \Pk\ 2 i n this term tells us that the stimulated 
creation process is enhanced by the same amplification factor as the spontaneous creation. 
Thus, the presence of initial quanta at early times will leave their imprint on the spectrum 
of perturbations at late times. This is important to take into account in computing the spec- 
trum of primordial perturbations resulting from inflation, as we will show in the remaining 
sections of this paper. 

It is worth noting that renormalization in the late-time Minkowski space does not affect 
the second term in square brackets in (13.91) . which is proportional to Tr[/x/VV]. So the potential 



impact of renormalization on the two-point function and the power spectrum [29H31|. and on 
higher order correlation functions do not interfere with the effects of the stimulated creation 
of perturbations. Therefore, we will not consider renormalization in the rest of this paper. 



IV. SCALAR POWER SPECTRUM AND STIMULATED CREATION OF 
PERTURBATIONS DURING INFLATION 

We are now ready to determine the observable effects on the CMB if there were pertur- 
bation quanta already present at a time to about 60 e-foldings before the end of slow-roll 
inflation. As we showed in the previous section, these quanta would give rise to the stim- 
ulated creation of quanta. As we have seen, in the absence of this stimulated creation of 
quanta, there would be no such observable effects because of the huge dilution of the original 
quanta during the 60 inflationary e-foldings. As we showed in the last paragraph of section 
ITT1 [see (12. 7p ], the presence of such quanta at time to is consistent with the scenario in which 
the unperturbed part of the inflaton causes inflation through its potential energy density 
function V(<p). We loosely refer to to as the initial time and the state at time to as the initial 
state. 

Inflation provides an elegant mechanism to generate inhomogeneities in the early 
universe 0. These inhomogeneities are believed to have originated as a consequence of the 
rapidly expanding inflationary geometry. This expansion induces the creation of quantum 
fluctuations of the inflaton field and also those of metric perturbations. These scalar and 
tensorial fluctuations evolve with time while their physical wavelength (in momentum 
space) is smaller than the Hubble radius Rh = H~ l during inflation. When their physical 
wavelength becomes larger than Rh (this is referred to as the exit from the Hubble radius), 
these fluctuations stop evolving and remain constant. When they re-enter the Hubble radius 
at a later stage of radiation or matter domination, they constitute the initial conditions for 
the inhomogeneities in the universe that will further evolve to produce the temperature flue- 
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tuations in the CMB and the LSS. We analyze here the production of scalar perturbations 
during inflation. The extension to tensor perturbations is straightforward. We consider 
the general case of single field slow-roll inflation. In this scenario, the homogenous part of 
the inflaton field, (f)(t), slowly rolls down its potential V(<p) towards a minimum during the 
evolution. The slow-roll parameters are defined as e = (Mp/2)(V' /V) 2 and rj = Mp(V" /V), 
where the prime here means derivative with respect to 4>{t) (Mp = 1/ V8ttG is the reduced 
Planck mass). In the slow-roll regime, defined as e,\i]\ <C 1, the evolution of the inflaton 
field is such that the potential V(<f)(t)) changes very slowly and takes the role of an 
effective cosmological constant in the Einstein equations. The universe then expands in a 
quasi-exponential way with H = a{t)/a{t) « ^(8^(^/3)^(0) remaining nearly constant, 
where a = d t a. The slow-roll parameters then can be approximated by e ~ l/(2Mp)(<p/H) 2 
and r) « -4>/(4>H) + l/(2Aff,)07# 2 ,. 

Scalar perturbations are usually characterized by the gauge-invariant quantity 7Z (the 
comoving curvature perturbation) which, during inflation, is given by 

TJ 

1Z = V + -^5<p , (4.1) 

where 5<j)(x,t) represents the inflaton perturbation field and ^ is a scalar perturbation of 
the metric, related to the curvature of the spatial metric through = 4V 2 \l//a 2 (t) (for 
details see, for instance, jsl, [HI]). To fix the time and spatial reparametrizations we choose 
the so-called comoving gauge in which 

5<P = , 9ij = a 2 (t) [e" 27 % + hij] , dih i:i = ht = . (4.2) 

In this gauge, the inflaton field is unperturbed and all scalar degrees of freedom are encoded 
in 71. In this section, it is sufficient for us to work to linear order in the perturbation 1Z. 
We expand the field in Fourier modes 

n(x, t) = J2 Mt)^ , #s(«) = ^n k (t) + A\_ s nut) . (4.3) 
k 

At linear order, the field lZ(x,t) is described by the action So = 1/2 f d 3 xdrz 2 [R, 12 — (cflZ) 2 ], 



where 1Z' = d T lZ, and the field modes obey the equation [32 



d 2 K k 2dzdK k 2 

~nr + ~i — — + kn k = o 1 (4.4) 

dr z z dr dr 

where r = J dt/a(t) is the conformal time and z = a<f)/H. In the slow-roll approximation 
z~ x dzjdr w (1 + 3e — tj)/t. The (normalized) solutions to this equation (considering e 
and 77 as constants) obeying the adiabatic condition [l[ (and the de Sitter symmetry for H 
constant) are 

K k {r) = (-Trr/ALh^H^i-rk) , (4.5) 

where Hjp is a Hankel function with index [i = 3/2 + 3e — 77. The amplitude of these modes 
goes to a constant value a few e-folds after the time r k at which the mode k exits from the 
Hubble radius (k/a(r k ) = H). It is this fact that makes the variable TZ interesting in the 
study of scalar perturbations. The spectrum of primordial scalar perturbations generated 
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during inflation is usually characterized by the so-called power spectrum P n . This quantity 
can be defined in terms of the momentum space two-point function evaluated a few e-folds 
after the Hubble exit time tj,. When the state of the field is given by the vacuum state 
[defined with respect to the modes (14.51) ]. we have 



<0|^Al°> = % 1+ 4,o^)> (4-6) 



where 5 



where the superscript indicates that the previous quantities are evaluated using the vacuum 
state. 

Let us suppose now that the state of the field is described by a statistical density operator 
p of the type defined in section [Til The power spectrum P^(k) in this case is defined by 

nP^A]^h 1+ ^,oPn(k) , (4.8) 

where a straightforward computation gives the following expression for P-j^(k). 

P&$) = Pnim + Tr\pN u ] + Tr[pi\L fe -]) . (4.9) 

Consequently 

Al(ky = Al(k)°(l + Tr[pN n ] + Tr\pN_$) , (4.10) 

where Tr \pN%] is the average number of quanta in the mode k present in the initial state 
and the superscript p indicates that a quantity is evaluated using the statistical density 
operator p. 6 We can observe in the above expression how the stimulated creation of quanta 
enhances the magnitude of the power spectrum by a factor (1 + Tr[pN^\ + Tr[piV and 
introduces the anisotropics that TrfpiVg] could contain. We consider now the constraints on 
the quantity TrfpiV^] coming from the present CMB observational data [tJ. 

First, we note that the high degree of isotropy of the CMB Q constrains the variation 
of Tr[/x/Vg] with the direction of k to be small. Therefore, it is a good first approximation 

to suppose that there is no dependence on the direction of k in the initial state so that 
Tr[pN^] « Tr[pN_^\. [However, if the evidence for anisotropies were to become stronger, 



5 In the limit L — >• oo this definition is consistent with the fact that the quantity A^(fc) is the contribution 
per cilnfc to the two-point function at coincident points, (0\lZ 2 (x, t)\0) = J dk/k A^(fc) . 

6 A similar analysis holds for tensor perturbations, and the result is P^ t (k) = P°(fc)(l + Tr[/9t7V|] + 
Tr[/9(7V* -]), where P^(k) is the power spectrum for tensor perturbations when the initial state of 
those perturbations is the vacuum state, pt is the statistical density operator describing the initial 
state of tensor perturbations at the time to, and Nt is the operator for the number of quanta of 
tensor perturbations in mode with coordinate momentum k. As a consequence, the tensor-to-scalar 
ratio takes the form r = Pfc* / P^(k) = r° ■ a, where r° = Ph/P^(k) = 16e and the extra factor 
a = (1 + Tr[p t Ni] + Tr[p t A^-])/(l + Tr[pN^] + Tr[piV_ £ ]) depends on the initial states of tensor and 
scalar perturbations. 
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then it would not be unreasonable to investigate possible anisotropies of the initial state as 
a potential source for anisotropies and other related anomalies (see, for instance 33| and 
references therein) that may be present in the CMB and LSS)]. In the rest of this paper, we 
consider only isotropic states of the initial perturbations. 

Next, we consider the constraints that come from the measured parameters A(k p ) and n s 
that characterize the observed power spectrum. The observed power spectrum is written in 
terms of these parameters as follows, 

Al(k) = A(k p ) , (4.11) 

where k p is a pivot (or reference) scale. The observations of 7-year WMAP (together with 
BAO and Hq Mean) [?J provide values for the two parameters, 

A(k p ) = (2.441±°;°^)1(T 9 , n s = 0.963 ± 0.012 , (4.12) 

where k p = 0.002 Mpc -1 . 

Our theoretical predictions now follow from (I4.10p . Using (14.71) and (I4.10p . one finds that 



and 

g i _ ^ h, (1 + 2 1Y W _ 

dink dink 
We now impose the observational constraints that 

A p (k p ) = A(kp) = (2.441+°;°g)10" 9 , n p s = n s = 0.963 ± 0.012 , (4.15) 

with the measured values in (I4.12p . 

Because the value of H 2 /e is uncertain at the time Tk p when the mode at the pivot scale 
exited the Hubble radius, the first condition in ( I4.15P does not strongly constrain the value of 
Tr[pNj: ]. On the other hand, the second condition in (14.151) imposes more severe restrictions 
on the dependence of TrfpA^] on k. In particular, this condition constrains the quantity 

dln< ~ 1 ^J^ pjVfc ^ to be not much bigger than the slow-roll parameters e ~ \rj\ ~ 1/100. 



V. PRIMORDIAL NON-GAUSSIANITIES FROM STIMULATED CREATION 
OF PERTURBATION DURING INFLATION 

As we have seen, the power spectrum alone has limited potential in revealing detailed 
information about the dynamics and interactions of the field or fields driving inflation. 
Higher-order correlation functions, such as the three-point function, are more sensitive to 
those aspects of inflation (see for instance (sj). In the following, we compute the three-point 
function in momentum space, to leading order in the perturbations, when the initial state of 
scalar perturbations is given by a density operator p of the type introduced in section [Til We 
obtain the values of the primordial non-gaussianites, those that were present after inflation 
ended, before their reentry into the Hubble horizon of the radiation- or matter-dominated 
universe. 
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The full computation of the three-point function using the vacuum state was given by 
Maldacena in |34|. Here, we generalize those computations to the case of a more general state 
described by a statistical density operator of the type discussed in section [TTJ As explained 
there, this density operator allows for the possibility that there were perturbations already 
present at a time to, about 60 e- foldings before the end of inflation. Our calculation includes 
the effects of the stimulated creation of perturbation quanta. This stimulated creation is the 
reason why the presence of such quanta early in the inflationary epoch cannot be ignored. 
We limit our calculation to the case of single-field inflation. 

Let us consider the three point function in momentum space for comoving curvature 
perturbation TZ(x), parametrized in the following way 



Tr[p% i (r)^ 2 (r)7t A ; 3 (r)] 



(5.1) 



where the Kronecker delta is a consequence of translational invariance and enforces the 
condition that ^ki = 0. The quantity Bn(ki, k 2 , k 3 ) is known as the bispectrum. As in the 
previous section, we work in the comoving gauge (14. 2p . In order to compute the bispectrum 
at leading order in the perturbations, we expand the inflationary action around the spatially 
homogenous solution to third order 



S = S (4 > ,g flu ) + S 2 ('Jl 2 ) + S 3 ('R. 3 



(5.2) 



The background action So(4>, g^), which depends on the background inflaton field 0(t) and 
the background metric, defines the Hubble rate H and the slow-roll parameters. The action 
S 2 = 1/2 J d 3 xdr z 2 \R! 2 — (dTt) 2 ], where TV = d T 7Z, defines the free evolution of the mode 
functions TZk(r). The action S3 defines the interaction Hamiltonian. The form of this action 
was computed in (34[. We refer the reader to [H| for further details. Maldacena obtained a 
simple form of the cubic terms in the action by using the following field redefinition 



K = K r 



+ 



2<pH SM 2 P H 2 



lit 



AMp H 2 



d- 2 (n c d 2 n c ) . 



(5.3) 



At the lowest order in the perturbations, TZ C is the same as TZ, so at that order the mode 
functions of 1Z C and 1Z are also the same. At leading order in the slow-roll parameters, 
S 3 = - J drH int , where 



<p x a 3 (r) ( -| 1 HM^ 2 n?d~ 2 Kl. . 



(5.4) 



We compute the three-point function at leading order in curvature perturbations and in 
the slow-roll parameters by using time-dependent perturbation theory. In the interaction 
picture, the three-point function in momentum space is given by 



Tr[p% i (r)^ 2 (r)^ 3 (r)] 



Tr 



P 



(5.5) 



int) ' 
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where T and T indicate the time and anti-time ordered product respectively, the superscript 
/ indicates operators in the interaction picture, and To is the conformal time corresponding 
to the cosmic proper time t at about 60 e-foldings before the end of inflation, as discussed 
earlier. Because comoving curvature perturbations Ttl remain constant outside the horizon, 
we can evaluate the previous quantity at late-times, when r — > 0. Additionally, in the 
computation we assume that ki ^ 0. The two terms in the last line of the previous equation 
can be computed explicitly (see Appendix for details) and the result is 

(0)^(0)^(0)]= (5.6) 

F(p, ki, k 2 ) + cyclic permut , 



<f>H 



+ 



4Mp H 2 



^3 



TO 



dr' Tr [p[TZl (0)7% (0)7% (0), H(Jf)\] = 

2 2 k\k\ 



(5.7) 



\^ l)fi P n(ki)P n{k2) 



M 2 H 2 kl 



G(p, ki, k 2 , k 3 ) + cyclic permut 



where 



F(p, h, k 2 ) = TV[p(2A^ + 1)(2A^ + 1)] - 6^ 2 Tr[pN Ri (N Ri + 1)] , (5.8) 
and G(p, k±, k 2 , ks) is a symmetric function under cyclic permutations of (k\, k 2 , k 3 ) given by 

(5.9) 



G(p, ki, k 2 , k 3 ) = 

(1 - cos (k t r )) 



kt 



Tr\pN^ N^ ] + TrfpiV^J + cyclic permut 



ki,k 2 



Tr [piV^iV^ + 1)] + cyclic permut 



+ 



COS (fciTp)) 



Trfp^^J + TtlpN^] - Tr[pN R N h ] + Tr[pA^] 



TrfpA^A^ + 1)] ) + cyclic permut 



where k t = ki + k 2 + k 3 and k\ = k t — 2ki. As we mentioned at the end of section HVJ we 
have assumed that Tr[pNz] « TrfpiV in the previous computations. (If it were necessary, 
it would be easy to introduce the effects of anisotropics.) In the previous expressions, the 
quantities H, <fi and <fi appear in combinations that can be expressed in terms of the slow-roll 
parameters, which can be regarded as nearly constant during inflation. Putting both terms 
together and taking into account (14.71) and (14. 9p . we obtain 



3e-2r) + e 



x 

.2 
'2 



Tr[p^(0)^ 2 (0)^ 3 (0)] = 5 (Efc - ))0 P£(h)P£(k 2 

"1 

2 

where we have defined 

F(p, ki, k 2 



(5.10) 



k 2 + k 2 \ - - k 2 k 2 - - - 

1 , 2 2 f(p, ki, k 2 ) + 4e-^ g(p, ki, k 2 , k 3 ) 



cyclic permut 



f(p, ki,k 2 ) 



Tr[p(2i\^ + l)]Tr[p(2iVg 2 + 1)] 



(5.11) 
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t T T T\- G(p,k 1 ,k 2 ,k 3 ) 

k >- h) = T r|P (2iv 5i + mM2N h + 1)] ■ (5 - 12) 

The denominators in / and g arise from the definition of P-^(k). Expression (I5.10p constitutes 
our result for the bispectrum when the initial state is given by the density operator p. 

If there are no particles (i.e., quanta of the perturbation field) present at time To, then 
Tr[piVV] = for all k, in which case we have / = 1 and g = (1 — cos {k t To))/k t . We can 
choose To such that at that time all the modes that we observe in the CMB are well inside 
the Hubble radius, so fcjr = ki/(a(r )H) ^> 1 and hence \k t T \ ^> 1. Averaging over the 
large range of uncertainty in this value, it is reasonable for us to take 1 — cos (k t r ) ~ 1. 



Then (I5.10p agrees with the result obtained in 34(, in which the pure vacuum state was 
considered. 

The presence of particles in the initial state introduces new interesting features. The 
most remarkable new element is the presence of the denominators proportional to ki present 
in the function G(p, k±, k 2 , k 3 ). The fact that ki can approach arbitrarily close to zero in 
some configurations opens the possibility of having a large contribution coming from those 
configurations, producing an enhancement of the non-gaussianities. This fact was already 
noticed for a pure initial state given by a Bogoliubov transformation of the Bunch-Davies 



vacuum m 



18[ [19[. There, the momentum configurations corresponding to flattened triangles 
for which two momenta are collinear were considered. However, as we will show, the most 
pronounced enhancement appears for configuration in which one of the three momenta is 
much smaller than the other two, the so called squeezed configuration. In the following, 
we particularize expression (I5.10p for each of the flattened and squeezed configurations, 
respectively. 

• Flattened configuration 

Let us consider first the case of flattened configurations, in which the momenta forming 
the triangle are such that k 2 ~ k 3 « —ki/2. For this configuration we have k t « 4/c 2 , k\ ~ 
and k 2 ~ k 3 ~ 2k 2 . 

In the case with no particles in the initial state, we have f(p,ki,kj^i) = 1 and 
g(p, ki, k 2 , ^3) = (1 — cos (k t r ))/k t « 1/kt, where, as before, we have averaged the cosine to 
zero. Equation f 1 5 . 1 j) reduces then to 

(01^(0)^(0)^(0)10) « 5 (E £ i) P^kjP^h) (31c - IO17) • (5.13) 

Let us consider now the case in which the initial state is described by the density oper- 
ator p. We can estimate the values of / and g by assuming that there are no significant 
correlations between initial particles with different momentum, so we can approximate 7 
Tr[pN%Ng] ~ Tr[pN^]Tr[pNg]. Additionally, if we assume that the average number of ini- 
tial particles is equal to or greater than one, and as we showed earlier, that the value of 
Tr[pNr] does not vary significantly with k, it follows that / is the ratio of quantities of order 
(Tr[pi\y) 2 and that f(p,h,k^i) ~ 0(1). 



7 For the density operator of the form given in (|2.1|) . this assumption implies that p(ni(ki), ?i2(fo), •■ 
p(ni(ki))p(ri2(k2))... . 
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On the other hand, because in this configuration fci ~ 0, the quantity g is essentially 
given by the term proportional to 1/ki, which has the form 

, r r r, 1 - cos&To) ^Ip^M ~ l^Wl) + Tr\pN Sl ] + lTr[pN %2 ] 
g{p,k 1 ,k 2 ,h) » - 



~h Tr[p(2A^ + l)]Tr[p(2N %2 + 1)] 

Pt 2 
ft l T 



2!ifci 



x 0(1) , (5.14) 



where we have expanded (1 — cos (&1T0)) for |&iTo| — > and have again used the fact that the 
factor involving the average number of initial particles is of order 1. Note the extra factor 
k\rQ introduced by the expansion of the cosine. Expression (15.101) then reduces to 

1^^(0)^(0)^(0)] » PUh)P£(k 2 )[0(e, V ) + 18e(fc 2 T )(fciT )] . (5.15) 

If we write the magnitude of the three-point function in the flattened configuration in 
terms of a parameter /jvz, in the form 

(^(0)^(0)7^(0)) = S^^PnfyPnih) f NL , (5.16) 

we obtain from (15. 13j) and (1 5 . 1 5 j) the following ratio of the value f^ L for an initial state 
described by the density operator p, to the value f^ L for the vacuum initial state 

|p ~ q , 18£ in (fc 2 r )(fc 1 r ) + 0(1) ~ (fc 2 r )(fc 1 r ) + 0(1) . (5.17) 

JNL 6le iU? 7 

Since k 2 is well inside the horizon at the early time To, we have that \k 2 To\ ^> 1. On the 
other hand, we have that k\To — > for the flattened configuration. We can conclude that 
there is not necessarily a large enhancement in the primordial non-gaussianities for flattened 
configurations when quanta are present in the initial state for the canonical interaction 



lagrangian. As shown in jl8j, different conclusions could be reached if higher-derivative 



interactions were considered. 
• Squeezed configuration 

Let us consider now the configuration in which two of the momenta forming the triangle 
are much bigger than the third one, for instance k\ ~ k 2 ^> k 3 . Then, we have k t m 2ki, 
k\ w k 2 « ks, ks « 2k\. The three-point function in the squeezed configuration is often 



parametrized in the following way |2l|,|22j 



(7^(0)7^(0)7^(0)) = 5 {E ^ Pn{h)Pn{h) l ^f NL ■ (5.18) 

In the case with no particles in the initial state, we have f(p,ki,kj^i) = 1 and 
g(p, ki, k 2 , fc 3 ) = (1 — cos (k t To))/k t — l/k t . Then, expression (I5.10p reduces to 

(01^(0)^(0)^(0)10) « 5 (E ^ i))0 PUh)PUh) (6e-2p) = 5 (E g i)i0 P^P^h) (1-nJJ) , 

(5.19) 
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where the scalar spectral index in the vacuum case is n = 1 — 6e + 2rj. Therefore, we have 
f% L = ^(1 — n s ), in agreement with previous analysis [34| . 

In the case where the initial state is described by the density operator p, we still have 
that f(p, ki, kj^i) ~ 1, provided the initial average number of quanta in the modes hi and kj 
is equal to or greater than one, and correlations between particles with different momentum 
are not large. Additionally, we have 

- r T > 1 - 22! <^ T ») 2Tr[piV gi W g J + Tr[ pJ V s J + Tr\pN h ] l^crafeTj) 
»(**.,*.,*») » s TM2N ^ ~ 1)mp{2Nh + 1}] ~ S xO(D ■ 

It is important to notice that, because k^ corresponds to a mode whose effects we wish to 
measure in the present CMB or LSS, its wavelength must have been well inside the Hubble 
radius at the early time r . That implies that k^TQ = k%/ (a{To)H) ^> 1, where H is the 
inflationary Hubble rate. Therefore, we can not expand the cos (A^To) above, and we have 

Tr[^ t (0)£ & (0)^(0)] » 5 (E ^ ))0 P^(fc 1 )P^ 3 )4e|(l - cos(fc 3 T )) , (5.20) 



from which we can obtain 



/ ^^S e | (1 " cos(fc3To))w r|' (5 - 21) 

where we have again averaged the cosine of large argument to zero. Note that j p NL is 
proportional to the slow-roll parameter e, so the size of the non-gaussianities is exactly 
zero for an exact de Sitter inflation. We can now compute, for the squeezed configuration, 
the ratio f^i/fNi between the magnitude of the non-gaussianities when there are particles 
present in the initial state and their magnitude in the case of the initial vacuum state 



f P NL w k i 4e ^ h 
Inl h 6e - 2r] k 3 



NL - (5.22) 



Because in this squeezed configuration we have k\ ^> k 3 , it follows that f p NL j ' f% L S> 1 • 

We can estimate the value of the ratio in (15.221) by considering the range of scales that 
we can measure in the CMB. The range of angular multipoles I for which we have some 
confidence about the measurements of the temperature fluctuations of the CMB (i.e., for 
which uncertainities coming from cosmic variance, Sunyaev-Zeldovich effect, etc., can be 
neglected) is approximately t G (10,2000). On the other hand, for £ ^> 1 the contribution 
of comoving momenta k to the multipole I is peaked arround k ~ £/ri, where Tl is the 
radial coordinate of the surface of last scattering. 8 Therefore, the maximum enhancement is 
achieved when ki w 2000/r L and k 3 ~ 10/r L . In that case ~ 200. It is important to note 

Jnl 

that this enhancement takes place when the number of initial quanta is equal to or greater 
than one, but does not depend on the exact value of the average number of initial quanta. 
This estimate is within the range permitted by the effects of backreaction, as discussed after 

(ETD. 



This is just a consequence of the properties of the spherical Bcsscl functions jeikr^) appearing in the 
relation between k- Fourier space and angular- Fourier space. For further details see, for instance, 0). 
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The analysis of this section shows that, when perturbation quanta are present in the 
initial state at the onset of inflation (or at an early time when perturbation wavelengths we 
can observe today were within the inflationary Hubble radius), the size of non-gaussianties 
produced during inflation would be considerably enhanced. The strongest enhancement 
would occur when the momenta appearing in the bispectrum form a squeezed configuration, 
for which there could be an enhancement (as compared to the initial vacuum state) of order 
one hundred in the parameter f^L that is usually employed to characterize non-gaussianities. 

It is important to mention here that our result contrasts with those reached in some 
previous analyses of non-gaussianities. In j23[ , they proposed a consistency relation for the 
three-point function of curvature perturbations in momentum space in the squeezed limit. 
Their consistency relation asserts that the size of the bispectrum in the squeezed limit is 
proportional to (1 — n s ), with the only assumption being that the inflaton is the only dynam- 
ical field. It follows from their consistency relation that, for those single- field inflationary 
models, the size of the parameter f^L that characterizes the primordial non-gaussianities 
in the squeezed limit must be of order of the slow-roll parameters. Any detection of non- 
gaussianities for the squeezed configuration would then rule-out any inflationary scenario 



with a single dynamical field. Recently, some weak aspects of the argument used in [23 
were pointed out in [35f . where a more rigorous proof of the consistency relation was at- 
tempted. However, in [35] they pointed out that they were not able to derive the consistency 
relation in all generality. They did succeed in proving it more rigorously for the vacuum 
initial state, but not for the general class of states we have considered in the present paper. 
Our results show that when particles are present in the initial state, there are new contribu- 
tions to the three-point function, or bispectrum, that come from perturbative interactions 
among the quanta produced by stimulated creation resulting from the presence of initial 
particles. These contributions can be as much as two orders of magnitude larger than the 
contribution coming from the vacuum. 



VI. CONCLUSIONS 

In this paper, we have analyzed the primordial spectrum of scalar perturbations produced 
during cosmic inflation when the initial quantum state is not necessarily the vacuum state. 
We have restricted our analysis to the scenario of single-field slow-roll inflation with the 
canonical Lagrangian. We have considered an initial state described by a statistical density 
operator p. In this mixed state, there are quanta of the perturbation field (i.e., particles) 
already present at the time t when the longest wavelength modes that we can observe today 
exited the Hubble radius during inflation. If this time is around 60 e-folds before the end of 
inflation, we do not need to consider trans-planckian excitations in p. The density operator 
p is formed from a mixture of pure states having definite numbers of perturbation quanta 
present at to- This density operator is natural in models of inflation in which our observable 
universe arises from the inflation of a small patch of a larger universe. 

Theoretical constraints can be imposed on the maximum average numbers of initial 
quanta, Tr [piVg], that are present in p for the occupied modes k. These constraints can 
be obtained by imposing renormalizability of the theory and demanding that backreaction 
effects from the initial quanta do not modify the background inflationary dynamics. These 
conditions restrict the average number of initial quanta but, as shown at the end of section 
\U\ do not necessarily restrict Tr[pA^] to be less than one. 

When p contains a non-zero average initial number of quanta, the phenomenon of stim- 
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ulated creation of perturbations occurs during inflation. We have computed the expression 
for the power spectrum and the bispectrum of scalar perturbations by taking into account 
the stimulated creation process. We have shown that the quanta present in the initial state 
do leave an imprint in the power spectrum and bispectrum of perturbations. Hence, present 
and future observations of the CMB and LSS can be used to constrain the initial state. In 
section HVl we showed how the expression for the power spectrum is modified by the presence 
of initial quanta. We used the observational data for the amplitude of the power spectrum 
and the spectral index to constrain the function Tr[pN^\. The isotropy of the CMB and 
the value of the spectral index severely restrict the variation of Tr[piVr] with k. However, 
due to our uncertainty in the value of the e slow-roll parameter and the Hubble rate during 
inflation appearig in the expression for the amplitude of the power spectrum, no significant 
constraints can be obtained on the magnitude of TrfpAM from the power spectrum. 

In contrast, the bispectrum (the three-point function in momentum space) is much more 
sensitive to the size of Tr[piVr] and constitutes a more suitable tool to probe the initial 
state. We have obtained, in section [V] the complete expression for the bispectrum when 
the initial state is described by p [expression ( I5.10p ], which includes the vacuum state as 
a special case. The most remarkable effect of the initial quanta is an enhancement of the 
bispectrum for certain configurations of the three different momenta involved. We have 
found that the largest enhancement appears for configurations in which two of the momenta 
are much bigger than the third one (for instance, fci ~ fc 2 > the so-called squeezed 
configuration. In this limit, for the range of perturbations accessible to observations, the 
size of the parameter f NL parametrizing the non-gaussianities can be as much as two orders 
of magnitude larger than its value for the pure vacuum initial state. 

This result is important for interpreting the observations made by the PLANCK satellite 
and by surveys of large-scale structure to determine the size of non-gaussianities. Our 
results indicate that, if non-gaussianities are found to be observable in the squeezed limit, 
then single-field inflation is not necessarily ruled out, as would follow from [23|. Instead, 
such an observation could be interpreted as a consequence of a non-vacuum initial state, 
and valuable information about this state could then be obtained from those observations. 
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Appendix: Explicit computations of the three-point function in momentum space 

In this appendix we present some details of the computation of the three-point function 
used in section [V] Let us consider that the initial state of scalar perturbations is described 
by the density operator p introduced in section [TTJ Additionally, as we mentioned at the 
end of section IIV} the isotropy of the CMB constrains the anisotropies contained in p, 
so we assumed that Tr[pAg] rs Tr[ pN g|. The generalization to the non-isotropic case is 
straightforward. Following Maldacena j34(, we employ the following field redefinition 



K = K r 



1 



2 0# 8MJ,# 2 
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Using this field redefinition, the interaction hamiltonian takes the simple form 

H int (r) = -Jd 3 x a(r) 3 (j^j HM p 2 K?d- 2 K' C . (6.2) 

In order to compute the three-point function in momentum space, we use time dependent 
perturbation theory at leading order in perturbations. Then we have 

^\pn %i {r)n h {T)n %s {r)] = (6.3) 



Tr 



Tr^(r)^(r)^(r)]-< / dr'Tr^^^r)^^)^^),^^)]] + 0{H, 



int) ' 



where T and T indicate the time and anti-time ordered product respectively, the superscript 
/ indicates operators in the interaction picture, and To is a conformal time about 60 e-foldings 
before the end of inflation. Because comoving curvature perturbations "RL remain constant 
outside the horizon, we can evaluate the previous quantity when r — > 0. Additionally, in 
the computation we assume that fc, 7^ 0. We will compute the two terms in the last line of 
the previous equation separately. 

The first term can be rewritten, by using the field redefinition (16. ip and using the convo- 
lution theorem, as follows 

1^^(0)^(0)^(0)] = 1^^(0)^(0)^(0)] + (6.4) 
1^^(0)^(0)^0)^0)1(51 -V) 2 + cyclic permut . 



1 

+ 



AMI H 2 B 

6 v 



The first term in the rhs of the previous equation vanishes, because of the form of p [see the 
discussion after ( 12. 61) ] . The trace appearing in the second term can be written, using (14. 3p . 

as 

^Tr[p^ gi (0)^ 2 (0)^(0)^ g3 _.(0)] =5>[p(^ fel (0) + A^T^O)) x (6.5) 
p p 

x (a^ 2 (o) + AK k ni(o))(A/iz p (o) + Al p ^;(o))(^ 3 _^ |& _ p1 (o) + ^ 7^(0))] . 



There are six (non-vanishing) expectation values in the previous expression that we have to 
compute. We show two of the calculations, as examples: 

E^^^V- ( fe-fl]^( )^( )^(°) w ft-si( )= ( 6 - 6 ) 
p 

2 ^\p A kAkAA) ~ ha^kAkAA)) l^(o)l 2 |7^ 2 (o)|\ v 



2Tr[p(% + 1){N %2 + 1)] - % life Tr[p%(A^ + 1)] 1^(0)1^(0)1%^,, 



+fc2+fc3,0 
2 I 



,0 ' 
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and 

^Tr[pA t _.A^A^ (g3 _ p) ]^ i (0)^ 2 (0)^;(0)^ |fe „ p1 (0)= (6.7) 

v 

= TripN^ + l)]|^ fcl (0)| 2 |^ 2 (0)|\, i+A : 2+A r 3i0 . 

By computing the four remaining expectation values, one obtains the result for expression 
(E2D as 

^Tr[p^ i (0)^^(0)7t c y0)^ 4 ^(0)] = 2F^k 1 M^kA^^ k M\ 2 h 1+ kM , 
p 

(6.8) 

where 

F(p, h, k 2 ) = Tr[p(2iV gi + 1)(2N %2 + 1)] - S^TclpN^N^ + 1)] . (6.9) 
In a similar way, we obtain the value of the last term on the rhs of expression ( 16.41) 



\ E ^KfeW^CO)^)^^)] A - = ( 6 - 10 ) 



Introducing (16. 8p and (16 . 1 0[) in (16 .4p . we obtain 

Tr[p^(0)^ 2 (0)^ 3 (0)]= (6.1i; 



^ 6 V0; n,(2fc 



It remains to compute the second term in the last line of (16. 3p . It contains the expectation 
value of the commutator of three fields operators and the interaction Hamiltonian. Let 
us compute the first term of the commutator. By Fourier expanding the field operators 
appearing in the interaction Hamiltonian and using the convolution theorem twice, we have, 
at leading order in perturbations 



— i 



jf dr'TV[^i(0)^ 2 (0)^ 3 (0)£r/ nt (r')] = -ijf dr'a\r') UA HM? L 3 x (6.12) 

x E E i^%(o)^ 2 (o)^^ • 

The previous expressions contains 20 expectation values similar to (16. 6p and (16.71) . but now 
involving six creation and annihilation operators. By computing those expectation values 
and by taking into account that, at leading order in slow-roll parameters 

n k (o)d T ,ni(T') = (S- ) V' , (6.13) 
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( I6.12p reduces to 



—i 



cir'Tr[p^ i (0)^ 2 (0)^ 3 (0)i// nt (r / )] = 8 {mfi 



H 4 {HY M 



4L 6 



k\k 2 k 3 \k\ k\ k\ 



I _ e i{ki+k 2 +k 3 )T 

ki + k 2 + k 3 



IV[p(^ i + l)(^ 2 + l)(iV^ + l)] 



(6.14) 



5t 



kito 



Tr [P%(% + + X )] + c y clic Permut 

i _ gi(-fci+fc2+fc3)n> 5r r 

[TtW^fe + ^ + !)] - + 1)]] 



+ 



-fcl + fc 2 + &3 



1 _ g*(- fc l- fc 2+fc3)T0 (5^ 

M^^Ws, + 1)] " [piV^ + l)(JV fe + 1)]] 



+cyclic permut I + 



-ki - k 2 + k 3 



I _ gi(— fci— fe2— ^3)7 

-A;i -k 2 -k 3 



+ cyclic permut ] + 



Tr^A^J - ( ^TrfpA^A^ + l)JVjJ + cyclic permut ) 



The other term of the commutator can be computed in a similar way producing 

r-O 

rlr'TrlnTT* <-r'\T^ . 

ki ^ ' fc2 kz 



dr' Tr [pHUr')nl (0)7% (0)&i 



S 



H A ( H\ 2 M~ 2 ( 1 

+ 



TO 



(Efc,),0 4L 6 \ ^ ) kl f, 2ks y k 2 ' k 2 ' k 2 



I _ e i(ki+k2+k-j,)T 

ki + k 2 + k 3 



Tr[pNrN u Nr]- 



fcl &2 ^3 
'5 



fcl,fe2 



Tr[pAr gi (Ar- + + cyclic permut 



k:'J 



+ 



J _ gi(- fc l+' c 2+fc3)ro 

-kx + k 2 + k 3 



[Tr[p(i\^ + l)JV fe JVj 



5 



&2>&3 



^[^ + 1)^(^ + 1)]] 



+cyclic permut J + 

1 _ gi(-fel-fc2+fc3)ro r 

_ fci _ fc2 + fcs [Tr^ + 1)(N R2 + l)JV fc ] - -^Tr[p(A^ + l)%iV^]] + 

+cyclic permut^ + 



J _ gi(— fcl— fe— fe3)T0 

-fci -k 2 -k 3 



Tr[p(% + 1)(^ 2 + 1)^3 + 1)] " 



-^Tr \pN %i (N ni + l)(iV^ + 1)] + cyclic permut J 



(6.15) 
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By subtracting (16.151) from (I6.14p . we have 



— ? 



dr> Tr [p p% (0)7^ (o)^ 3 (0) , hUt')]] 



(6.16) 



5 



( Efc»),o 
L 6 



HV H A 



i>j 



G(p, fci, k 2 , k 



3) , 



where G(p, fci, A; 2 , A^) is a symmetric function under cyclic permutations of (kx, k 2 , k 3 ) given 
by 



G(p, kx, k 2 , h) = 

(1 - COS (k t T )) 



(6.17) 



TrlpN^N^] + TrfpiV^J + cyclic permutj + 1 — 



Tr[pN^ (N^ + 1)] + cyclic permut J 



+ 



+ 



1 — cos 



>s(fciT )) / 



npN^} + Tr[pN % Nj: 3 ) - Tr[pN^N %3 ] + Tr[^J 



Tr[pN^ 2 (N^ 2 + 1)] ) + cyclic permut 



with kt = k\ + k 2 + k 3 and ki = k t — 2k{. 

Putting (16. lip and (I6.16P together and taking into account (14. 7p and (14. 9p . we obtain 



^7^(0)7^(0)7^(0)] = 5^ P£(fci)P£(fc 2 ) 



x 



(6.18) 



- (^3e - 277 + e * fc2 2 j /(p, fci, fez) + 4e-^#(p, k x , k 2 , k 3 ) 



cyclic permut 



where we have defined 



f(p, h,k 2 ) = 



F(p, k ly k 2 



g(p, h,k 2 ,k 3 ) = 



Tt[p(2N %1 + l)]Tr[p(2iV^ + 1)] 
G(p, h, k 2 , k 3 ) 



Tr[p(2iV gi + l)]Tr[p(2iVg 2 + 1)] ' 
Expression (I6.18P is the three-point function that we presented in section IVl 
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